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For simulation studies of (macro) molecular liquids it would be of significant interest to be able 
to adjust or increase the level of resolution within one region of space, while allowing for the free 
exchange of molecules between open regions of different resolution or representation. We generalize 
the adaptive resolution idea and suggest an interpretation in terms of an effective generalized grand 
canonical approach. The method is applied to liquid water at ambient conditions. PACS numbers: 
02.70.Ns, 61.20.Ja, 80.82, 05.10.-a, 05.20.Gg 
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Many (macro-) molecular systems display phenomena 
and properties that are inherently multiscale in nature, 
and thus particularly challenging for both theoretical and 
experimental investigations. Here computer simulations 
represent a powerful tool of investigation, though truly 
large all-atom simulations of complex molecular systems 
are neither possible, because of CPU requirements, nor 
in many cases very useful, as they often produce an ex- 
cess of data. To the aim of linking large scale generic 
properties and local specific chemistry, a variety of scale 
bridging simulation techniques, ranging from hierarchi- 
cal parameterizations of models covering different levels 
of resolution, see e.g. [IkZI, to interfaced layers of differ- 
ent resolutions (see e.g. |8l4l8|) have been developed. 

As most techniques are sequential - the whole system 
is treated on one level of resolution at a time - switch- 
ing between resolutions may be performed only for the 
whole system. In many cases however it would be desir- 
able to adjust the level of resolution on the fly only in 
a smaller well defined region of space, i.e. considering 
more details, according to the problem of interest, while 
keeping the larger surrounding on a coarser, computa- 
tionally much more efficient level. This idea has been 
successfully employed for problems of crack propagation 
in hard matter [l9[. However for soft matter or liquids 
inherent molecular number fluctuations pose special chal- 
lenges. They require a simulation setup, which allows 
molecules or parts of large molecules to cross bound- 
aries between areas of different representation without 
any barrier, keeping the overall thermodynamic equilib- 
rium intact. While crossing between regions of differ- 
ent resolution, they either become more coarse-grained 
or resume more atomistic detail that is they acquire or 
lose degrees of freedom (DOFs). In this context, adap- 
tive resolution MD simulations would represent a natural 
Grand Canonical set up: An atomistic region/cavity is 
coupled to a large, ideally infinite coarse-grained region 
acting as a particle/molecule reservoir. In general for 
such a coupling scheme the only requirement is that the 



coarse grained surrounding supplies/removes molecules 
and heat in a way that in the atomistic region (U,V and 
T remain constant. The adaptive resolution method 
(AdRcsS, Adaptive Resolution Scheme), which allows 
one to couple two systems with different resolution within 
one MD simulation describes a first step in this 

direction. Molecules change their resolution and thus 
their number of DOFs when moving through space and 
passing from one region through a transition zone into 
the other region. The method has successfully been ex- 
tended into the quantum [2^ and the continuum 2^ 27 1 
region. A different approach based on the interpolation 
of Lagrangians can be found in Ref. [2^ In the follow- 
ing we will generalize the AdRcsS concept and show that 
the atomistic/ fine grained subsystems can be viewed as 
a system in contact with a (huge) coarse-grained par- 
ticle/molecule reservoir, thus effectively representing a 
/iVT ensemble on the level of the high resolution sub- 
system, while the overall systems remains to be an NVT 
ensemble. Actually, by linking the coarse grained region 
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in addition to a continuum, as shown in Refs. [2; 
this is easily extended to a ^VT ensemble on the level of 
the whole system. We now first shortly review the adap- 
tive resolution method AdResS and the concept of ther- 
modynamic force which assures overall thermodynamic 
equilibrium. We then extend this to a more general cou- 
pling scheme and apply it to the example of an adaptive 
resolution simulation of liquid water, where the coarse- 
grained water model matches the compressibility of the 
all- atom model, while the pressure is higher by three or- 
ders of magnitude. 



In the AdResS method molecules smoothly change 
their level of resolution as a function of position by mov- 
ing through a transition region as illustrated in FiglT] 
In order to accomplish the process described above, two 
different representations are interpolated by a transi- 
tion function w{x) with w = 1 in region A, the fine 
grained more microscopic or all-atom one, and w ^ Q 
in the coarse-grained region B, respectively. Further 



2 




FIG. 1: Pictorial representation of the adaptive simulation 
box and local molecular representation. On the left, is the 
low resolution (coarse-grained) region, indicated by B, and on 
the right, is the high resolution (atomistic) region A. In the 
middle is the transition (hybrid) region H with the switching 
function w{x) (curve in grey), 

let US assume that a particle in B is the coarse-grained 
representation of a collection of atoms in region A and 
the position of the coarse particle is the center of mass 
of the corresponding atomistic molecule. With this 
in mind two molecules a, /? in region A (fine-grained, 
w{Xa) = ^{Xp) = 1, Xa, Xjj being the x-coordinates 
of the centers of mass of the two molecules, cf Fig. [T]) 
interact via the standard forces between all individual 
atoms of the molecules; in region B (coarse-grained, 
w{Xa) = w{Xfi) = 0) two coarser molecules interact 
with each other via central forces; just as in a bulk sys- 
tem of all atom molecules or coarse grained molecules, 
respectively. In all other cases of mixed resolution in- 
teractions, these two forces are interpolated as illustrated 
in Fig.[T]so that the total force between two different 
molecules consist of two parts. The all-atom contribution 
reads: 

Ftp = w{X^)w{Xp) ^ a^p (1) 

While the coarse-grained parts is: 

J-f^ (1 - w{X^)w{Xp)) Ff^ a^p (2) 

V^ p, is the force between atom i of molecule a and 

atom j of molecule j3 and F^^ is the coarse-grained cen- 
ter of mass - center of mass interaction. Note that in 
the transition region the intra-molecular atom-atom in- 
teractions are left unaffected. For technical details of 
the implementation we refer to the supporting mate- 
rial. Altogether we simulate a system at constant vol- 
ume, partitioned in two different (ideally macroscopic) 
regions (coarse-grained and fine-grained) by a transition 
region of finite width. This region facilitates the exchange 
of molecules and ensures equilibrium between the coarse 
and fine grained regions, moreover, by construction (as 



explained below), should not affect the density fluctua- 
tions (compressibility) . Beyond the crucial technical task 
of assuring the thermodynamically correct coupling be- 
tween the two large regions, the transition region does 
not have, a priori, any thermodynamic/physical meaning. 
An important characteristic of the adaptive approach is 
that while the overall number of molecules is conserved, 
this is not the case for the number of atoms in the in- 
dividual regions. Eventually we aim for a flat density 
profile throughout the whole simulation box, including 
the transition region, so that we have a thermodynamic 
state point well defined without (kinetic) barriers for the 
exchange of molecules. 

In the same spirit, the temperature has to be well de- 
fined and equal throughout the system. However, while 
this is simple to achieve in the "bulk" regions, the tran- 
sition region requires some care since in this region a 
change in the number of degrees of freedom occurs. By 
switching from coarse-grained to atomistic or vice versa, 
degrees of freedom (DOFs) are reintroduced or removed 
and the related equilibration is accomplished by a " local" 
thermostating procedure (for details see the supporting 
material) . The approach presented above, paves the way 
to an open system simulation scheme within a generalized 
Grand Canonical framework. 

A common method to derive coarse grained interac- 
tions is based on matching the center of mass radial 
distribution functions between the molecules in the two 
systems. As a result the compressibilities k remain 
unchanged poj by construction while the pressures in 
the fine and coarse-grained representation can be differ- 
ent [131 . In terms of the grand potential —pV this reads 
(for identical volumes) 

PA{tJ'A,T)V ^PB{jJ.B,T)V; ka^kb, (3) 

where /i^ and /is are the respective chemical poten- 
tials. This will result in density variations in the sys- 
tem. Instead we would like to have a constant density po 
throughout the whole simulation box, while keeping iden- 
tical compressibilities. Different pressures at the bound- 
aries of the transition region will produce a drift force 
on the molecules, leading to an unphysical inhomoge- 
neous system. Different compressibilities would affect 
the molecule number fiuctuations in the - still finite - 
regions. Thus we keep ka = kb and propose to exactly 
compensate this pressure generated drift force by a ther- 
modynamic force Fthix) acting on the molecules and Eq. 
[3] can then be rewritten as (see Fig.[T]) 

(^PA{t^A,T) + ^ Fth(x)dx^ 1/ = Ps(Mb,T)I/ 

(4) 

with pq being the reference, uniform, molecular density. 
EqH] tells that the subsystem A (B) is in equilibrium 
with a reservoir ps {pa) respectively, notwithstanding 



3 



the fact that the pressures, pa and Pb, and local chemi- 
cal potentials, and /x^, rnay be different. To overcome 
these variations, the related thermodynamic work, which 
exactly compensates any related drift force between the 
different regions, is provided/absorbed by Fth- Let the 
all-atom system, at a given density po, serve as refer- 
ence system, then the force on a molecule with mass Ma 
balancing —Vp{x), reads 



Fth(a:) 



Ma 



Vp(x) 



(5) 



where p{x) is the local pressure at the target density 
Po- Fth is applied for a < x < h, where a (b) is the 
left (right) boundary of the transition region. In EqlHl 
p{x = a) = PA and p{x = b) = pb- Now the density will 
remain unchanged. With this generalization the original 
equation for the coarse-grained forces in the transition 
region, Eq. [2] is changed so that the force acting on 
the center of mass of molecule a reads: 



The above line of arguments does not require the pres- 
sures in the two ("macroscopic") regions, pA and ps, 
to be the same, allowing to couple rather different sys- 
tems. Formally this simply means that the thermody- 
namic force as well as the thermostat can perform work 
on the molecules to adjust the virial pressure and the 
thermal energy, respectively, while passing from one re- 
gion to the other. This procedure allows one to couple 
almost arbitrary systems and keep them in equilibrium 
with respect to each other. For typical coarse-graining 
procedures like (iterative) Boltzmann inversion or Re- 
verse Monte Carlo which are mostly based on radial dis- 
tribution functions and reproduce structural rather than 
thermodynamic aspects very well, this is of special ad- 
vantage. 

Unfortunately the optimized pressure profile p{x) is not 
directly accessible as it has to be measured under the 
constraint of enforcing the equilibrium density po ui the 
system. Thus we developed an iterative approach based 
on the density profile. Within a linear approximation a 
first estimate of the local pressure p(x^p{x)) for an en- 
forced overall constant density po is: 



p{x,p{x))) KipA 



1 



PqKt 



{po - p{x)) 



(7) 



where kt is the (constant) isothermal compressibility, 
and p{x) is the, non uniform, density to which the system 
would adjust itself to without applying an external force. 
In order to enforce the uniform density po , Eq. [7] and 
Eq. [5] allow to obtain an initial thermodynamic force, 
that can be refined iterativcly through 



Fl+\x)=Fl^{x) 



Mg 



Vp'(x), 



(8) 



until the system evolves to the target uniform density. 
Here for i = 1, we have p^{x) = p{x) and Fl^{x) = 0. 
The prefactor 1/pqKt can be interpreted as a variation 
of local chemical potential, by the identity [2^ 



dp 
'dp 



1 



V,T 



P KT 
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a fact explored already before in this context [23| . In the 
following we will demonstrate the power of this concept 
for the illustrative case of liquid water, where the com- 
pressibility of the SPC/E water and a structure based 
coarse-grained model are the same, while the pressure in 
the coarse-grained system is 6000 barfsi}. 

Coarse graining liquid water with a focus on the liquid 
structure usually leads to models for which the pressure 
is significantly higher while the compressibility remains 
essentially unchanged. Thus pressure correction terms 
are often used, which add linear terms to the coarse- 
grained potential. The compressibility is then affected as 
(6) Kj- cx / r^{g{r) — l)dr and small deviations at high r val- 



ues can have dramatic effects [31|. However, for adaptive 
resolution methods, but also solvent-solute systems, it is 
important that the compressibility in the coarse-grained 
and all-atom regions is the same. This minimizes finite 
size effects for our (in most cases on purpose as small 
as possible) atomistic region, imposed by the - still fi- 
nite - CG region. For many applications in computa- 
tional biology, but also for studies of water itself, the 
SPC/E model is employed [111- A standard structure 
based coarse-grained model perfectly fits the radial dis- 
tribution functions, but does not comply with the tetra- 
hedral packing 3l|, |33[ . It has -within the error bars- the 
same compressibility but a very high pressure of about 
6000 bar, making it an ideal test case for the thermody- 
namic force principle described above. If not properly ac- 
counted for, this pressure difference between the different 
models would lead to density distortions throughout the 
simulation box and create a transient flow of molecules 
from the coarse-grained to the atomistic region. Using 
at the outset of the iteration scheme of Eq. [5] as 
described before, after only four iterations for the ther- 
modynamic force wc arrive at a flat density profile, as 
illustrated in Fig. [2l 

The thermodynamic force creates a stationary state 
with the correct target density, removing both the effect 
of the pressure difference of the models on the distribu- 
tion of molecules in either region and the artifacts from 
the change of resolution in the transition region. The 
advantages are also clearly displayed, when looking at 
local density fiuctuations. While in the thermodynamic 
limit the molecule number fluctuations are related to the 
compressibility through pUbTut = l-iv)^^ ' here 
analyze an effective compressibility in terms of molecule 
number fluctuations in slabs along the axis of resolution 
change. The increased compressibility for a pressure cor- 
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FIG. 2: (a) Density profile for the water simulation for the un- 
corrected AdResS simulation of water and several iterations of 
the thermodynamic force with A being the all-atom SPC/E 
region and B the coarse-grained region. The final density 
profile shows deviation of less than 0.4 %. (The arrow indi- 
cates the development of the density with the iteration.) (b) 
Compensation of the local pressure change between regions 
by the thermodynamic force, (c) Molecule number fluctua- 
tions in slabs (width Aa; = Inm) along the x-axis for pressure 
corrected and uncorrected coarse-grained models. 



rected coarse-grained SPC/E water model leads to an in- 
crease in the local density fluctuations. Within our setup 
we do not alter the effective compressibility anywhere in 
the system as shown in Fig. Actually the molecule 
number fluctuations in subvolumes of region A, are the 
same as in similar subvolumes of a full all-atom simula- 
tion. Hence the particle distribution in region A shows 
the expected Gaussian distribution as it is the case in the 
full all-atom simulation, shown in Fig. [31 Leveling out 
any pressure undulations in the transition region, while 
keeping the compressibilities the same also means that 
there is no barrier for molecule exchange between regions 
(see supporting information) . Coupling two systems with 
different pressures means that the thermodynamic force 
performs work on the molecules, when they cross from 
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FIG. 3: Particle number distribution in a slab of width / = 
4 nm and for an atomistic zone in adaptive resolution with the 
same width. Both distributions match the expected Gaussian 
behavior with mean /x = {N) and variance — phsTKriN) 



one region to the other. In practice the integral over the 
thermodynamic force over the transition region does not 
have to be zero, which would be the case if the pressures 
on either side were the same. Thus for the present setup 
the thermodynamic force has to exactly compensate the 
local pressure variation due to the constant density. The 
comparison of the final local pressure (see supporting in- 
formation for details) and the thermodynamic force as 
derived in form of an external local pressure field is shown 
in Fig. [2b- They perfectly match, illustrating the over- 
all consistency of the approach. Thus by applying the 
thermodynamic force we arrive at a robust and efficient 
simulation setup, which enables us within a single MD 
simulation to keep rather different systems with quite 
different intrinsic pressures but identical compressibili- 
ties in equilibrium with each other. 

The fact that the number of molecules at a given res- 
olution fluctuates suggests to view each subsystem as a 
reservoir of molecules for the other. Strictly speaking, 
to interpret such a setup as a Grand Canonical ensem- 
ble would require macroscopic reservoirs. In our case it 
is obvious that the number of molecules, and thus the 
reservoir, is finite. In fact, we have shown that the av- 
erage and the fluctuating number of molecules in each 
subsystem and even in smaller subregions throughout 
the whole simulation box are the same as analyzed in 
a full atomistic system. For the general case where the 
all-atom as well as the coarse-grained system are at the 
same (or even well defined different) densities the sim- 
ulation setup leads to a situation, where a difference in 
the excess chemical potential fi (the temperature part 
for the internal degrees of freedom is taken care of by 
the thermostat) between the regions, is compensated by 
the work performed by the thermodynamic force. Thus 
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on the level of the (large) subsystems we run a (/x, V, T) 
MD simulation, while we have a {N, V, T) ensemble glob- 
ally. The extreme example of liquid SPC/E water at 
temperature 300K and pressure 1 bar, where the coarse- 
grained model at the same temperature and density has 
a pressure of more than « 6000 bar, shows the power of 
the present approach. In addition, as mentioned before 
26l 27 1 one can link the coarse-grained region quite easily 
to a continuum domain, making it globally Grand Canon- 
ical. This opens completely new pathways to study sys- 
tems, where a rather high level of detail is only required 
within a rather restricted area, while at the same time 
the exchange with the surrounding is essential. First ex- 
amples of applications are (small) organic molecules such 
as amino acids or peptides in mixed solvents, e.g. water 
and urea. For peptides the conformation is expected to 
strongly couple to the urea density close to the peptide. 
There the present setup facilitates very fast and efficient 
high precision calculations of Kirkwood Buff integrals to 
estimate solvation free energies. More generally, among 
many others, potential applications include growth and 
structure formation phenomena in soft matter where lo- 
cally the assembly process needs details of the molecular 
structure (templated crystal growth, bio-mineralization 
etc.) and/or a controlled flux of additional molecules is 
provided. In conclusion we have derived a thermody- 
namically sound, technically robust, and easy to imple- 
ment way to have, in a single MD simulation one (small) 
system coupled with a different (large) system in Grand 
Canonical equilibrium. We exemplified this for the case 
of an atomistic system in Grand Canonical equilibrium 
with a source of atomistic molecules provided by a bath 
of coarse grained molecules. This is possible due to the 
work provided by the thermodynamic force in the tran- 
sition zone, acting like an active membrane, facilitating 
the transformation of particles/molecules. 
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I. FORCE CALCULATION 

In Eqn. 1 and Eqn. 2 of the original paper the interac- 
tions between molecules a and /3 subject to the interpo- 
lation of forces is written in a rather general form as 



for the center of mass force and 



w{X^)w{Xp)Y,¥ 

i,3 



A 



(1) 



(2) 



for the atom level interactions. From there one arrives at 
the interpolation formula for the total force between two 
molecules used in previous publications: 

= w{X^)w{Xp)¥tp + (1 - w{X^)w{Xp))Y'^f, (3) 

From this it is possible to see that, (i) the force F^^ 
in the hybrid region is actually a function of all the 
atomic coordinates {ri)i^a and {rj}j(zp and hence need 
not to be in the direction of the connection vector Rap 
of the centers of mass and can lead to e.g. rotation 
of the molecule, (ii) Since the intra molecular forces 
do not contribute to the total force on a molecule, we 
do not have to include the intramolecular interactions 
in the interpolation scheme; they are kept acting on 
the individual atoms throughout the atomistic and 
transition regions just as they are. 

In general Eqns. 1,2 describe four different cases for 
the forces between two different molecules, which apply 
to the whole simulation box: 

(i) For w(Xq,) = w{Xp) = two coarse-grained molecules 
interact via the effective coarse-grained interaction 



(4) 



leading only to a translational motion of the molecules 
which are now reduced to point particles. In this regime 
there is no intermolecular or intramolecular atomistic 
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contribution. 

(ii) For ■w{Xa) = w{Xj3) = 1 both molecules are in the 
atomistic region. Thus atom i of molecule a interacts 
with all atoms j of molecule /? 



./3 



J6/3 



A 



(5) 



The situation is that of a standard all-atom simulation 
with forces acting on each atom. 

(iii) For < w{Xa) < 1 and w{Xp) 7^ at least molecule 
a is in the transition regime and /? is not in the coarse- 
grained region. Then both parts of the interpolation con- 
tribute to the motion of the molecule, hence the coarse- 
grained part is also contributing to the force between the 
atoms. Following Eqns. 1,2 the force on atom origi- 
nating from all atoms of molecule /3 reads: 



il-w{XMX(^))j^F%. 

(6) 



The first part describes the standard atom atom inter- 
actions, while the second contains the distribution of 
the coarse grained interaction onto the individual atoms. 
By this the coarse grained interactions are naturally in- 
corporated into the all atom resolution, which we keep 
throughout the whole transition region, 
(iv) For w{Xa) and w{Xp) — only the second part 
of Eqn. [6l the center of mass contribution, survives and 
has to be distributed to the atoms of molecule a, which 
gives 



(7) 



for atom a^. 

With this scheme all forces between atoms and molecules 
are well defined throughout the whole simulation box. 
As mentioned in the main text, the thermostat needs 
some attention. Kinetic energy and the related tem- 
perature are both well defined in the coarse-grained and 
the atomistic region. More attention is required in the 
transition region. By switching from coarse-grained to 
atomistic or vice versa, degrees of freedom (DOFs) (in- 
cluding their related thermal energies) are introduced or 
lost. The related contribution to the thermal energy will 
be furnished by a standard local thermostat as has been 
discussed in Refs. (7l-[l0}. Thus the molecules/particles 
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are "handed over" properly thermalized at the bound- 
aries of the transition region. In summary, one can ap- 
ply any standard local thermostats like, e.g, a Langcvin 
thermostat, which automatically take care of the "latent 
heat" [1, 111 related to the fading DOFs. AdResS simu- 
lations from the point of view of the whole system cor- 
respond to simulations of a NVT ensemble, where N 
counts the number of molecules, but not the number of 
DOFs (which in our case are not strictly but only in av- 
erage proportional). In the current implementation we 
use a Langcvin thermostat, which acts without any in- 
terpolation on the atoms. The friction is proportional 
to the masses of the atoms and thus by themostatting 
each atom, the center of mass of every molecule obeys 
the fluctuation-dissipation-theorem as well [3]. A pure 
center of mass thermostat would also be sufficient, but 
slightly more difficult to implement in the integrator. 
The consequences for the definition of kinetic energy etc. 
that follow from "partial" degrees of freedom have been 
extensively discussed in 

The above discussion considers molecules which are re- 
placed by their centers of mass in the coarse-grained re- 
gion. However a more general mapping from explicit 
to coarse-grained, like a atomistic polymer mapped to a 
chain of spheres (so called "coarse-grained" beads) con- 
nected by springs, is possible. Then the interpolation is 
applied to the coarse-grained beads. 



II. PRESSURE CALCULATION 

We can numerically analyze the local pressure in the 
stationary state and its behavior when an external field 
is applied making use of the method proposed by Todd, 
Evans and Daivis The time averaged local pressure 
tensor is given by 
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^ (v„ •e^)(v„ •e^)M„\ (8) 

{q:|X„-2:|<Ax} ' ' 



where A is a cross-sectional area perpendicular to the 
x-direction. The first term corresponds to the ideal gas 
contribution calculated in a small slab centered at x with 
thickness 2Aa;, and the second represents the contribu- 
tion of the forces. Both sums run over molecules (which 
are labeled with a as before), where Mq,, Vq and Xq are 
mass, velocity and position. Consequently, the ideal gas 
contribution of the center of mass is fully included in the 
sum regardless of the weighting function of the molecule, 
since those degrees of freedom fully contribute to the av- 
erage in any representation. The internal degrees of free- 
dom, however, are not considered here. Therefore, the 
resolution of each molecule is implicit in the force con- 
tribution that depends on the weighting function. In the 
case of two-body forces, the second sum is reduced to 



the count of the forces that cross the area A. Particles 
far away from this plane (at distances bigger than the 
maximal cutoff of the non-bonded interactions) do not 
contribute to the sum. 



III. SIMULATION PROTOCOL 

The simulated system consisted of 8507 water 
molecules in a box of dimensions 16 nm x 4 nm x 4 nm at 
ambient condition, the simulation was performed with 
a modified GROMACS package Q]- Electrostatic in- 
teractions were treated with the reaction-field Q ap- 
proach using a cut-off length of 1.2 nm. At a temper- 
ature of 300K and a pressure of 1 bar we observe for 
a system of 8507 water molecules a compressibility of 
K«t = 4.56 X 10-i°toV^I3- 

The coarse-grained potential was obtained by mapping 
each water molecule in the center of mass to a single 
coarse-grained bead by an iterative Boltzmann inversion 
based on radial distribution functions. lOO steps of 0.2ns 
each were carried using the VOTC A software package 0] ■ 
As mentioned before, such a coarse-graining does not pre- 
serve the state point, yielding a coarse-grained pressure 
of approximately 6200 bar, while the compressibility is 
not altered. That is the local slope of the equation of 
state is conserved by such a procedure, while the abso- 
lute value is shifted. Additional corrections to reproduce 
the atomistic pressure can be applied 0, but have been 
shown to alter the compressibility [1] by a factor of 5, 
which one has to avoid in order to obtain the correct par- 
ticle number fluctuation in the different regions of space. 
To obtain the highest accuracy in the final density pro- 
file the thermodynamic force was extended by half the 
cut-off length (0.6 nm) into the atomistic zone (because 
of the mixed atomistic coarse grained interaction of two 
molecules, when one molecule is in the transition region). 
Without this extension the density deviated at the most 
by about 1 % from the target density, showing a drop in 
the atomistic zone. This is due to the fact that close to 
the coarse-grained zone, the atomistic molecules interact 
with the hybrid molecules an thus have a slightly altered 
state point. 



IV. DIFFUSION PROFILE 

To proof that molecules indeed diffuse back and forth 
between both regions without regard to the local level 
of resolution, we can measure diffusion profiles along the 
axis of resolution change. At t = to a-U particles in a slab 
of width Ax ~ Inm are selected which arc then "traced" 
(in terms of their density profile) over time (fig. [T]). This 
is done for a set of particles in both regions. As shown 
by the symmetry of the decay for both cases, the tran- 
sition region does not introduce any further transition 
barrier. On the coarse-grained side the particle distribu- 
tion decays significantly faster, illustrating the enhanced 
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diffusion constant in this region. 



8 9 10 11 12 13 14 15 16 



X [nm] 

FIG. 1: Interdiffusion between the two regions. This shows 
the density profiles over time for particles originally in a cubic 
subvolume with width Ax = Inm at t = io- The profiles are 
asymmetric as difi'usion constants are different in each repre- 
sentation. By this it is clear that both domains constantly 
exchange particles. 
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